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Abstract 

I propose a method for studying the product of bandwidths for the Harper-Hofstader 
model. This method requires knowledge of the moments of the midband energies. I 
conjectured a general formula for these moments. I computed the asymptotic repre- 
sentation for the product of bandwidths in the limit of a weak magnetic flux using 
Szego's theorem for Hankel matrices. I then give a first approximation for the edge of 
the butterfly spectrum and discuss its connection with P. Levy's formula for Brownian 
motion . 



1 Introduction 

A system of electrons on a square lattice in an uniform magnetic field displays an energy 
spectrum with multifractal properties. Many papers were devoted to this problem. See, 
for example, [§ ^ |g 0, 0, 0, ^, ^, |^, |2§. A lot of effort has been concentrated 



to find the Hausdorff dimension of the spectrum at a given magnetic flux. For a rational 
magnetic flux, i.e. $ = 27rp/g, with p and q relatively prime integers, there are q bands in 
the spectrum. The scaling hypothesis, [|12[, states that the majority of the bands behaves 



asymptotically as q~^, when q ^ oo, the rest of the bands being exponetially small. The 
Hausdorff dimension does not depend on the asymptotic behaviour of the exponentially 
small bands. The multifractal analysis is concentrated on the bands of the scaling type. The 
exponentially small bands were studied from a semiclassical aproach in a number of papers. 
Consult |]lT| and references therein. 
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The aim of this paper is to study the asymptotic properties of the product of bandwidths, 
so that both scahng and exponentially small bands will manifest themselves in the asymptotic 
formula. 

First I discuss a method to obtain the band spectrum. Second, I define more precisely 
the problem I will address. 

The Hamiltonian of electrons on a square lattice in a uniform magnetic field, in terms of 
magnetic translations Ti and T2, is given by 



where Ti and T2 obeys the following commutation relation 



i<I>A 



TiTs = e^'TsTi. 



In ( |1.2| ), $ is the magnetic flux. See [^, |2B[ for details about magnetic translations, 
the case A = 1 will be considered hereafter. 



For a rational magnetic flux. 



$ = 27r 



p 



the Hamiltonian becomes a q x q matrix. From now on, p and q will be relatively 
integers. 

One representation for the translation operators is 



where 
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Wi 
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W2 = diag{e^'^, e*^*, ■ ■ ■ , e" 
and 6j, j = 1,2 are real numbers (Bloch parameters). 

In this representation, the Hamiltonian is represented by the following matrix: 



H{e) 



( ai hi 
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where 



an{Q) = 2A cos($r2 + 62) {n = 1,2,- ■ ■ , q) 
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Chamber's relation, 0], states that the characteristic polynomial of the Hamiltonian 
matrix decomposes into a 9 -independent polynomial Pp/q{E) of degree q and a function 
h{9). 

det{H{e) -E) = P,/q{E) + h{e) , (1.8) 

h{9) = {-iy-'2{cos{q9i) + cos{q92)) . (1.9) 

When 6 varies, the function h{6) varies between —4 and 4, so, for a given rational flux 
$ = 27[p/q, the bands are obtained by intercepting the graph of the polynomial Pp/q{E) with 
the horizontal lines drawn at 4 and —4. Figure 1 shows the procedure for p = 1, q = 3. 
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Figure 1: Bands for P1/3 = —x^ + Qx 



For the case in the figure 1 there are three bands, because g = 3. Counting from left these 
are Bandi = [—1 — a/3, —2], Band2 = [1 — V^, —1 + V3], and Bands = [—1 + a/3, 1 + a/3]. 
Denote the widths of these bands Ai, A2,and A3 respectivelly. In general, for a polynomial 
Pp/q{E), the widths will be named, counting from negative to pozitive energies, Ai, ■ ■ ■ , Ag. 
The roots of the polynomial Pp/q{E) are eigenvalues for the Hamiltonian, for special values 
of Bloch momenta 9 = (vr/g, vr/g). Let us call these roots Ei,- ■ ■ ,Eg, counting from the 
negative to pozitive values of E. For the example in figure 1, p = 1, q = 3, these eigenvalues 
are: Ei = —a/6, E2 = 0, E3 = \/Q. The energies Ei, E2,- ■ ■ , Eg are called midband energies. 
These energies were studied in p|, ^ using Bethe-ansatz. 

The problem can be now stated: find the asymptotic formula for the products of the 
bandwidths, i.e. 

nA„. (1.10) 

n=l 

In the last section of the paper, I obtain that, for a week magnetic field $ = 271- -^ 0, 
the asymptotic formula can be arranged as a product of terms, each term approaching zero 
faster than the next neighbour term, (|5.14| ): 



HA. 



^ q ^exp(— gdo 



rfi + 0(l)). 
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Here, do and di are some constants. We will compare now this result, (|1.11|) , to results from 
multifractal analysis of the band spectrum. The multifractal analysis requires the power 
sums of the bandwidths to be evaluated and the goal is to find the function F{P) defined by 



,F{I3) 



T.K 



;i.i2) 



n=l 



where /3 is a real parameter. Taking the derivative with respect to /3 at /3 = in (|1.12|) we 
get 

F'(0)glng = ^lnA„ . (1.13) 



n=l 



Compare ( |1.13| ) wih the logarithm of ( |1 . 1 1| ) and you see that the term 2~'^ , which is due 
to exponentially small bands, does not appear from ( |1.13| ). To incorpoarate this term, we 
need to know how many bands out of the total number of bands ( which equals q) are 
exponentially small. Introduce a number / between and 1, such that fq is the number of 
bands that are exponentially small. Than the total number of bands can be decompose as 
q = fq+ {I — f)q and (1 — f)q is the number of bands that asymptotically behaves like q~^. 
Now the formula ( |1.12| ) can be rewritten as 



fqe^GW 



f)qQ 



Fm-i 
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n=l 



Here G{l3) is a function which describes the effect of all exponetially small bands. Than by 
taking the derivative with respect to /3 at /5 = 0, but now in ( |1.14| ), we get (using G(0) = 
and F(0) = 1. ) 



fG\0)q' + (1 - f)F\0)q\nq = ^ 1^ A 



9^ 

ra=l 
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In this case both leading terms in ( p. . 1 1| ) are present in ( |1.15D and by comparing the coeffi- 
cients we can write 



;i-/)F'(o) = -i , 

/G'(0) = -ln2 



;i.l6) 



The above formulae, ( |1.16| ), are valid for the case of a weak magnetic flux ($ — > 0). It will 
be interesting to flnd the fraction / as a function of the flux. By the scaling hypothesis, the 
fraction / is close to 0. I want to emphasise that the formula ( [1.14| ) is a hypothesis and work 
remains to be done to understand the behaviour of J2n=i A^. See |jl2| for a detailed study of 
multifractal properties and the bibliography therin. 

Before ending this section, I write the polynomial Pp/q{E) as a characteristic polynomial 
of a tridiagonal matrix. In this way it will become easier to compute the moments of the 
roots of the polynomial Pp/q{E), which we need to obtain the asymptotic formula ( |1.11| ). 
The moments will be the main object of study in Section 3. From this point to the end of 
the Introduction, the formulas are taken from |1^ . 



The roots oi Pp/q{E) are special eigenvalues and correspond for those 6 for which h{6) = 0. 
Studying the polynomial Pp/q{E) becomes more facile if we choose another representation 
for the Hamiltonian matrix : 



In this case the Hamiltonian matrix acquires the form: 



;i.i7) 
:i.i8) 



/ oi 
0,1 a2 



m 



02 



V a„ 





dn-i ; 
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with 



Chamber's relation in the above representation is: 

det{H{e) - E) = Pp/,{E) + h{e) , 

where 

h{9) = (-l)«-^2(cos(g^i) + (-If cos(g(^2 + ^i))) • 



;i.20) 
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This last representation is convenient because for the special value ^o = (^iOj^2o) = 
(0, 7r(l + p/q)), simultaneously 

^(^o) = (1.22) 

and 

o,(^o) = , (1.23) 

which makes H{6q) a tridiagonal matrix. Out of this tridiagonal matrix, recursive equations 
can be deduced for generating the polynomial Pp/q{E). Let Pq{E) , Pi{E) , ■ ■ ■ , Pq{E) be a 
sequence of polynomials generated by the following relations: 



where 



P^{E) = -E , 

Pn{E) = -EPn^i{E) - l3n-lPn-2{E) (u = 2, 3, ■ ■ ■ , g) 

/3„ = 4sinM— n) (n = 1, 2, ■ ■ ■ , g) , 



;i.24) 
;i.25) 
;i.26) 

;i.27) 



$ = 2'Kp/q . 



It follows that Pq{E) = Pp/q{E) and that the polynomial Pp/q{E) is the characteristic poly- 
nomial of the following tridiagonal matrix, which contains only positive numbers: 



M 



/O /5i 

1 [32 
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Pp,q{E)=det{M-E) . 
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2 Bandwidths' Product 



From the paper IJT^, it is known that the bandwidths obey the inequalities (with e = exp(l) 
2.718 ■■•): 



2(l + v^ ^ 8e 



n = 1 



) ; 'i ; 



(2.1) 



where Pi {En) is the derivative of Pp/q{E) computed at the eingenvalue En, which are the 
roots of Pp/q{E). As a consequence, 



(2 + 2v^)^ 



nUPLJE., 



< il An < (8e) 
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(2.2) 



The product 11^=1 Pp/qi^n) can be reexpressed as: 

n (E, - E,f . 

The formula ( p.3[ ) can be written as a product of two Vandermonde determinants 

1 



(2.3) 



n (E, - E,r 
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E, 
El 



■ Et^ 

1 E') Er) ■ ■ ■ E2 



1 El Ei 
^2 -E'2 



^9-1 ^r ... Ei~^ 

p/q 
S2k = JZEf, fc=l,2,3--- . 

n=l 

Putting it all together, it results: 



I Eq El ... E^' 



(2.4) 



By this observation, the product nn=i Pp/giEn) can be written in terms of the eingenvalue 
moments, defined as 

(2.5) 



(2 + 2^5)^^ < n A„ < 



:8e) 
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where cr„_i is the Hankel determinant: 



cr„-i 
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Sq+1 



Sq~l Sq Sg+i • • • S2q-2 



(2.7) 



To get the asymptotic representation of the product of the bandwidths, as g — > oo, a 
separate study of the eigenvalue moments, S2k, is necessary. The next section is dedicated 
to this study. 



3 Eigenvalues' Moments 



In the begining of this section, I am going to show that 



[fc2/4] 

S2k = q Yl '^i(^) cos{2njp/q) if q > k 

j=0 



(3.1) 



where aj are integer numbers that do not depend on p and q. Here [A;^/4] is the integer part 
of fcV4. 

To see this, note that the eigenvalue moments can be computed from the trace of the 
powers of the matrix M, ( [L.28 ): 

S2fc = Tr(M2^) . (3.2) 

Let us consider two examples, for fc = 4 and k = Q : 



g-i 



TriM^)= ^(2/?^ + 4/3„_i/5„ 



m=l 



y-i 



TriM'^) = J2 i^Pi + QPl-iPm + Q(3m-iPl + 6/3„_2/5^^-i/5„ 



(3.3) 



(3.4) 



m=l 



where /3m = if m < 0. The sum can be extended to q, because Pq = 0. We get then (here 
C := e^'^'P/''): 

<? 
S4 = Tr{M^) = J2 (2(2 - C"" - C-"")^ + 4(2 - C""-^ - C-™+^)(2 - C"' - C""')) (3.5) 

rra=l 

= 4g(7 + 2cos(27rp/g)) . 



In general. 



S2k — 22 Z2 di^^h-ljPm-j' " Pm~lP\ 

m=l lo+h^ \-lj=k 



m 1 



(3.6) 
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where j takes values between and k — 1; compare with ( p^) . The coefficients diQi^...i. do 
not depend on p and q. They only depend on k and on (Zq, h- ■ ■ Ij). Moreover, we can write 

/ . Pm-j ■ ■ ■ Pm-lPm ~ Z_^ 9n^ " ; (,'J- ' J 

m=l n 

where f^^ are integer numbers which do not depend on p and q and the index n runs from 
over a finite range. In the above sum, the only nonzero terms will come for those powers C,n 
which do not depend on m. This fact is only a consequence of the structure of (/q, ^i, ■ ■ ■ , ij) 
and not on p and q. The sum over m will collect q identical terms so in the result the 
parameter q will be factor out. Note that for a fixed k the same (/o,/i, ■ ■ -/jO-structure is 
present for all q > k. For q < k, some (/q, h- ■ ■ Ij) are not present. For example, for q = k 
the term Pm-k+i ■ ■ ■ Pm-iPm is not present, whereas for all g > fc is present. I end up here 
the motivation for ( |3.1| ). The maximum value for j in ( |3.1| ) is [/c^/4]. This will be explained 
later, in section 4, after formula 



To proceede further, note that in the formula for S2k, (|3rT|), the parameter p appears only 
in the ratio p/q. Because of this, the analysis of the eigenvalue moments can be extended 
to all real magnetic fiuxes. Indeed, for a irational a;, let us consider a sequence of rational 
approximants p„./g„ — ;> u;, as n — >■ cx). For each n we get 

SA,n = 4g„(7 + 2 cos(27rp„/g„)) . (3.8) 

So, in the limit n — > oo, 

< E\uo) >:= lim ^'=^ ^'■" =4(7 + 2 cos(27rcu)) , (3.9) 

where -Ej>, 3 = !■ ■ ■?« are the special eigenvalues (roots of Pp^/q^^{E)) for the fiux $„ = 
2TTpn/qn- Here are the first five limit eigenvalue averages: 



<E^{uj) > = 4 , (3.10) 

<E^(cj) > = 4(7 + 2 cos(27rcj)) , (3.11) 

<E6(w) > = 4(58 + 36cos(27rtu) + 6cos(47rtu)) , (3.12) 

< E%uj) > = 4(539 + 504 cos(27rcj) + 154 cos(47rcj) + 24 cos(67rc<;) + (3.13) 
4cos(87ra;)) , 

< E^^Lu) > = 4(5486 + 6580 cos(27rcj) + 2770 cos(47rcj) + 780 cos(67rcj) + (3.14) 
210 cos(87rc<;) + 40 cos(107rcj) + 10 cos(127rc<;)) . 



Now, a natural question conies to mind: is there a general formula for S2k, or equivalently 
for < E^f'iuj) >? 

At this time, I have a partial answer to this question. From the first eigenvalue moments 
the following statement can be infered: 



[fcV4] 



2 'tC' (2fc-2m-l)!! k\ 



< E^\-) >= ( 1^) 1 + E ^ ' r,:v^u\. \,,^ p-(^^ -^"^(--) ' (3.15) 



m=l 



{2k-l)\\ (k~3^-2)V 



where for each m the positive integer number j^ is that one for which : 

2[(j„ + l)V4] + 2<2m<2[(j„ + 2)V4] . (3.16) 

The first pairs (2m, j^) are (2, 0), (4, 1), (6, 2), (8, 2), (10, 3), (12, 3), (14, 4), (16, 4), (18, 4). 
Pm.{k) is a polynomial of degree 3m — jVn "" 3. This is so because the power of k in front of 
sin^"^(7rp/g) is 2m. There is one more detail about Pm{k) which can be found by inspection 
of the first few terms. Namely, the coefficients of the largest power in fc, i.e. \m in : 

P^(A;) = A„P™-J™-3 + . . . , (3.17) 

are generated by 

\/2t °° 

^ ^ = V \„,x^^ . (3.18) 

smh(v2x) „=o 

The following is a list of the first three polynomials Pm- The first eight polynomials are listed 
in the Appendix. 



Pi{k) = -\ , (3.19) 



P2{k) = ^{7k^-Ak-lh) , (3.20) 



P^{k) = (31A;^ - 60fc^ - 101A;2 + 66fc + 280) . (3.21) 



From the formula ( |3.15| ) we get S2k that we need in ( |2.7| ). Namely 

S2k = q< E^\p/q) > if q>k . (3.22) 

In order to compute the asymptotic representation of the Hankel determinant (p.7| ), a knowl- 
edge of Pm{k) will suffice, because the largest moment needed in (|2.7| ) is S2q-2 which is still 
given by the formula ( p.22| ), since g > (2g — 2)/2. Because there is no a general formula for 



the polynomials Pm{k), the asymptotic representation for the bandwidths' product for all 
magnetic fiuxes cannot be computed. 



Still, the case $ = 27ri when q —>■ oo can be analized, which is done in the last section. 



.9 



Before that, the eingevalues' moments can be further studied, due to summation formula 



( p. 17] ). The anlysis from the next section, besides a value of its own, it is also useful for 
understanting the asymptotic formula for the bandwidths' product for the case of a week 
magnetic flux. 



4 The Edge of the Butterfly Spectrum and Levy's For- 
mula for Brownian Motion 

In order to obtain a first approximation for the edge of the spectrum, we can make use of 
the general formula which gives the eigenvalue moments (|3.15| ). Let us write the formula for 
the eigenvalue moments as a power series in k and retain only the largest power in k. Call 
the result the first approximation to the eigenvalues' moments, and use (1) as a superscript 
to mark it: 



Q' 



,[kV^] 1 



E E^T = In ( T) E A^^A:- sin- vr^ ) . (4.1) 



For fc ^ oo we get 



E E]T = ?n(C,\)^— -^-^f^ , (4.2) 

j=i smh ( /esm(7r^ 



and 



-Pn' 



,.^ I , o ksm(TC- , . 

E^iu) = hm hm I g„(C2\)^— -y-^fV\^ 1 , (4.3) 

smh ( fcsm(7r^^) 



Pn/qn^'jJ k-*CO 



from which, the first approximation of the edge of the butterfiy spectrum is given by 

EiKc.) = 4e-^-(--) , (4.4) 

where < u; < 1 and the fiux is $ = 2Tiuj. 

Figure 2 presents, for each p/g, g = 53 , p = 1, 2 ■ ■ ■ 52, the spectrum Ej{p/q), j = 1 ■ ■ -q 
on a horizontal line. The vertical coordinate is called Flux ( though we represented only 
p/q and not 27rp/q). The horizontal coordinate is obviously called Energy. Superimposed 
on the above described spectrum is the graph of 4e~2^™('^'^) as a function of cu, which is 
the first approximation to the edge of the spectrum. Second order approximations can be 
obtained by summing the next largest power of k, namely all terms that contain fc^™"^. It 
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Figure 2: Spectrum for g = 53 and the function 4e 2'^™('^'^), the first approximation for the 
edge 



is quite clear that a full understanding of the sequence of polynomials Pmik) will enable us 
to decipher the fractal character of the edge of the spectrum. 

At this point it is worthy to mention works connected to the eigenvalues' moments for- 
mula, (|3.15|) . The geometric interpretation of the trace of the Hamiltonian can be found, for 
example, in the paper 0. For this, consider the square lattice Z^, where Z are the integer 
numbers. Let T2k be the set of all paths of lengths 2k which start and end at the origin 
(0, 0). For such a path F, let area(F) be the oriented area enclosed by F, see |^ for details. 
Then we have: 

Tr{H^^) = q Yl e'"^^'^'^^^) . (4.5) 

rerafc 

Compare this expression with the formulas for the eigenvalue moments. For example, for 
k = 2we obtain that 4-7 = 28 represents the number of paths of lenght 4 starting and ending 
at the origin, which enclose zero area. With the interpretation ([4.5|) for the eigenvalues' 
moments, it is easy to explain the maximum value [fc^/4] reached by j in ( p.l|) . The maximum 
value for the area(F), when F has a fixed perimeter of 2k, is achieved when F is a square 
with the edge of length k/2 (when k is even), so the area is k'^/4. By the same token, for k 
odd, the maximum area is [A;^/4]. 

As /c — > oo the random walks with 2k steps, aproaches (using a suitable renormaliza- 
tion, see 0), the 2-dimensional Brownian motion. The probability distribution of the areas 
enclosed by a planar Brownian motion was computed first by P. Levy, |^ . Levy 's result is 



E[exp{ig area)] 



g area/2 
sinh((yf area/2) 



(4.6) 



where -E[ ] is the expectation value and g is a real parameter. Also, in a recent paper ^1 
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Mingo and Nica studied the power sums of the areas. What the authors found is that: 

('')"' ^ (area(r))2™ = i?2™(fc), if k > 2m , (4.7) 



rerafc 



where R2m{k) is a rational function in k. The degree of R2m{k) (i.e., the difference of the 
degrees of the numerator and the denominator of R2m{k) ) is equal to 2m, and the leading 
coefficient of R2m{k), call it i'2m, is generated by: 



1^2 



^- = V-^(2z)2™ . (4.8) 

sin;z j^(2m)!^ ' ^ ' 

Both these results are connected with the generating function for the leading coefficients of 
the polynomials Pm{k), (|3.18| ). The formula for the eigenvalue moments, ( p.l5| ) go beyond 
the Levy's formula. From this perspective, a full understanding of the polynomials Pm{k) 
will put in a new light the planar Brownian motion. 



5 Asymptotic Representation of the Product of Band- 
widths for a Weak Magnetic Field and Szego's The- 
orem 

We now aim to compute the asymptotic representation of the product of the bandwidths, 
nn=i ^n, for a weak magnetic flux, $ = 27ri, when g — > oo. We will use the inequality 
( p.6|) together with the asymptotic representation of the Hankel determinant, (|2.7| ). The 
entries of the Hankel determinant are the moments S2k-, see (|3.22|) . Up to the flrst order of 
approximation, ( [4.2|) , the moments are: 

,^,s2 A;sin(-) 

S2k = q{'^) ■ , ,, -V.,, • (5-1) 

^ ^ smh(fcsm(-)) 

To go further, and because of lack of knowledge of the polynomials Pm{k), I assume that 
( for the case of a week flux), the second and all other orders of approximation are much 
smaller than the flrst order. Morover, because 

ksm(-) k 
Ci< . , ,, .'',,.. e^<C2, fc = l,2,---g-l , (5.2) 

smh(A;sm(-)) 

I will consider that the influence of the flrst approximation term, i.e. 

^'^ (5.3) 



sinh(A;sin(-)) 

on the asymptotic formula is of the same order of magnitude as given by the term e '^i . In 
( ^.21) Ci ~ 7re/sinh(7r) = 0.739 and C2 ~ e^/'^/sinh(l) = 1.169. Thus the moments S2k to be 
used to flnd the assymtotic formula for the Hankel determinant are: 

S2fc = g('fc')%-* , (5.4) 
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and it follows that 
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(5.5) 



The dimension of o"„_i is g. To find the asymptotic representation of crg_i, as g ^ cxd, 



we shall employ a theorem of G.Szego. One version of this theorem, |T^, gives the asymp- 
totic behaviour of a determinant Dn-i{f) '■= det(6^j,)o<;i,i/<n-i whose generic elemet can be 
expressed as 



-'l^u 



x"+V(x)rfx , 



(5.6) 



that is to say, the entries of the determinant can be expressed as the moments of a function 

fix). 

The theorem stays that, as n ^ oo. 



Dn^iif) = ^n-i(l) exp inco + g Eic| + 0(1) 
where c^, k = 0,1 ■ ■ ■ , are the Fourier coefficients of the function ln(/(cos0)) : 



(5.7) 



'' = 2^ 



-ike 



\n{f {cos 9))de 



{5.t 



To apply this theorem we need the function /. The generic term in the Hankel determinant 



aq-i can be expressed as an integral as followes: 



'J III/ — o 



X' 



>'+''K\x)± 



X 



where K {x) is the complete elhptic integral of the first kind: 



K(x) 



IT 



il-x''^ne)-'/'d9=^,F,i^,^,l;x' 
2 \2 2 



1 1 



Here x^ + x ^ = 1 and 2-^1 (ct, b; c; x) is the Gauss' hypergeometric series 

The function / is then 

fix)=K'ix) . 



(5.9) 



(5.10) 



(5.11) 



The determinat Dn-i{l) can be computed exactly and is given by the following product, 
using formula (2.2.15) in ||23[| : 



n-l 



Z}„_i(l) = 2n [k 



k=l 



_(^k 

2k^2k 



(5.12) 
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For large value of n 

Z}„_i(l) = 2-"("-i)+«(i) . 

Combining with the result of Szego's theorem, the asymptotic representation of crq_i 
reads as: 

a,_i = g^e-'^+i (^-^ j 4''('^-i)D,_i(1) exp gco + - ^ Jc? + ^(1) ' (5-13) 

Here Cq = 0.729. 

From this we get, as g -^ oo, for a weak magnetic field ($ = 27ri), 
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E 

ra=l 



n A„ ~ 2-'^\-'^ expi-qdo + rfi + 0(1)) , (5.14) 



where do, di are constant which cannot be made precise because of the innequality we started 
with, (|2.6|), and because of the approximations we used, (|5.1| ) and ([5.4|). 

Suppose now that the bandwidths are written as an exponential: 

An = e-f'"'^ . (5.15) 

From the above definition we obtain the asymptotic formula for the average value Hn '■ 

^-^^"=ln2 + ^-^ + 0fl') . (5.16) 

At the limit: 

<li>:= lim ^"=^^" =ln2 . (5.17) 

6 Conclusions 

This paper contains asymptotic results on the bandswidths in the Hofstadter spectrum for 
2-dimensional electrons in a magnetic field. Besides definite results, it also leaves some open 
questions. The important results are: asymptotic formula (^.14|) , the formula for the edge 



of the spectrum ( [4.4| ), and the conjecture on the general formula for the eigenvalue moments 
(S). 



The open questions are: What is the general formula for the polynomials P.in{k)1 How 
to obtain the fractal structure of the edge of the spectrum from the eigenvalue moments? 
How to find the asymptotic representation for the product of the bandwidths for every flux, 
using Szego formula? 
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8 Appendix 

I list here the first eight polinomials Pm{k), which appeared in eigenvalue moment formula 
( P.15| ). The following formulas are conjectured from the first 25 eigenvalue moments ( 3.101 ). 



Pi = -3"^ (8.1) 

P2 = 90"^(7A;^-4A;-15) (8.2) 

P3 = -1890"^(31A;^ - 60k^ - lOlP + 66k + 280) (8.3) 

P4 = 37800"^ (127fc^ - 1044A;^ + 2246/c^ + 328A;^ + 7k^ - (8.4) 

12244A;2 - 220A; + 25200) 

P5 = -3742200"^ (2555A;9 - 31887A;^ + 137946A;^ - 243774fc^ + (8.5) 

290499^:^ - 796527A;^ + 647416A;^ + 1661436A;2 - 
95184fc - 3991680) 

Pe = 10216206000-^(1414477A;^2 - 32155043fc^^ + 300761927A;^° - (8.6) 

1517115007A;^ + 4720753473fc^ - 10693488621^;^ + 20222071853A;^ - 
25923159133P + 11996284390fc^ + 11841173876A;^ + 
45697115160^;^ + 12156530400A; - 90810720000) 

P7 = -10216206000"^(286685A;^^ - 8595813fc" + 110951109/^^^ - (8.7) 
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819736201A;" + 3945073115A;i° - 13604875839P + 
36023850727A;^ - 73323185763A;^ + 108429092220A;^ - 
120643403168P + 122711359264A;^ - 65709283536A;^ - 
63479070720A;2 - 19852750080A; + 186810624000) 

Ps = 20841060240000^^(1 18518239A;^^ - 5242493627/^^*^ + 104055035258A;^^ - (; 
1233596237660A;^^ + 9844400576738A;i3 - 56735776354034A;i2 + 
247874192827336A;" - 847184326535620A;^° + 2294971824780007A;^ - 
4924658969019451A;^ + 8359347851375254A;^ - 11207926976651080A;^ + 
11409224332806816A;^ - 7994766402003888A;^ + 4489556615965152A;3 - 
3442931101232640A;=^ - 806985556876800A; + 4668397493760000) 
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